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Abstract In this paper we extend the theory of close encounters of a giant 
planet on a parabolic orbit with a central star developed in our previous work 
(Ivanov & Papaloizou 2004, 2007) to include the effects of tides induced on 
the central star. Stellar rotation and orbits with arbitrary inclination to the 
stellar rotation axis are considered. We obtain results both from an analytic 
treatment that incorporates first order corrections to normal mode frequencies 
arising from stellar rotation and numerical treatments that are in satisfactory 
agreement over the parameter space of interest. These results are applied to 
the initial phase of the tidal circularisation problem. We find that both tides 
induced in the star and planet can lead to a significant decrease of the orbital 
semi-major axis for orbits having periastron distances smaller than 5-6 stellar 
radii with tides in the star being much stronger for retrograde orbits compared 
to prograde orbits. Assuming that combined action of dynamic and quasi- 
static tides could lead to the total circularisation of orbits this corresponds to 
observed periods up to 4-5 days. 

We use the simple Skumanich law to characterise the rotational history of 
the star assuming that the star has its rotational period equal to one month 
at the age of 5Gyr. The strength of tidal interactions is characterised by circu- 
larisation time scale, t ev , which is defined as a typical time scale of evolution 
of the planet's semi-major axis due to tides. This is considered as a func- 
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tion of orbital period P bs, which the planet obtains after the process of tidal 
circularisation has been completed. 

We find that the ratio of the initial circularisation time scales corresponding 
to prograde and retrograde orbits, respectively, is of order 1.5-2 for a planet 
of one Jupiter mass having P b s ~ four days. The ratio grows with the mass 
of the planet, being of order five for a five Jupiter mass planet with the same 
P or b- Note, however, this result might change for more realistic stellar rotation 
histories. Thus, the effect of stellar rotation may provide a bias in the formation 
of planetary systems having planets on close orbits around their host stars, as a 
consequence of planet-planet scattering, which favours systems with retrograde 
orbits. The results reported in the paper may also be applied to the problem 
of tidal capture of stars in young stellar clusters. 

Keywords Tidal interactions • Dynamic tides • Rotating stars • Planet 
formation and evolution 

1 Introduction 

One explanation for the existence of closely orbiting extrasolar giant planets 
is through the mechanism of orbital migration, induced through the gravi- 
tational interaction with the protostellar disc (e.g. Lin & Papaloizou 1986). 
Notably this mechanism is expected to produce planets in circular orbits that 
are coplanar with the the stellar equator. However, the presence of high or- 
bital eccentricities amongst extra solar planets is indicative of a strong orbital 
relaxation or scattering process. This led Weidenschilling & Mazari (1996) and 
Rasio & Ford (1996) to consider planet-planet scattering of planets initially on 
nearby circular orbits that resulted in one planet being scattered onto a highly 
eccentric orbit that had a close approach to the central star. Subsequent tidal 
interaction with the central star was then postulated to lead to the forma- 
tion of a Hot Jupiter. Papaloizou & Tcrquem (2001) investigated the strong 
dynamical relaxation of a planetary system with 5 to 100 planets starting on 
closely packed circular orbits. This type of relaxation also readily produced 
planets in a highly eccentric orbits with close enough approaches to the cen- 
tral star for stellar tides to lead to orbital circularisation. The evolution of the 
planet orbit due to combined action of planet-planet scatterings and tides has 
also been considered (e.g. Nagasawa et al 2008, Marchi et al 2009). 

More recently measurements of the inclination of the orbital planes of tran- 
siting protoplanets through the Rossiter-McLaughlin effect has led to the dis- 
covery of a significant number of close orbiters on highly inclined or even retro- 
grade orbits (see eg. Winn et al. 2010, Hebrard 2011, and references therein). 
This indicates that a planet-planet interaction that can lead to high inclination 
orbits with high eccentricity that are subsequently circularised has been active 
in some cases at least. Thus tidal circularisation of orbits with arbitrary incli- 
nation with respect to the central star need to be considered. Tides induced in 
the planet are expected to be independent of orbital inclination because of the 
rapid attainment of pseudo synchronisation (see Ivanov & Papaloizou 2007). 
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However, this is not the case for tides induced on the star. For example in a 
frame rotating with the star, tidal forcing from a retrograde orbit will appear 
to have a higher forcing frequency than tidal forcing from a prograde orbit. Ac- 
cordingly the tidal response will differ. This leads to a larger transfer of energy 
from the planet orbit to the star during a single periastron flyby (see also Lai 
19970 This type of phenomenon may modulate the distribution of observed 
Hot Jupitcrs on inclined circular orbits that may have been circularised with 
tides acting on the central star playing a significant role. 

It is the purpose of this paper to study the tidal interactions of planets 
on highly eccentric orbits, with orbital planes arbitrarily inclined to the stel- 
lar equatorial plane, taking account of tides induced on the central star and 
estimating their contribution to the circularisation process. In making these 
studies we follow the approach of previous papers (Ivanov & Papaloizou 2004, 
2007; Papaloizou & Ivanov 2010) and approximate the initial stages of orbital 
circularisation as being through a sequence of encounters of the planet on a 
parabolic orbit with the star. We find typical time scales, t ev , of evolution 
of the planet's semi-major axis a due to dynamic tide exerted in the planet 
and in the star, as a function of the period P b s , which the planet has after 
the circularisation has been completed and the planet has settled on a low 
eccentricity orbit. This period can be obtained from observations. 

It is important to note that the aspects of dynamical tides considered 
in this paper can only efficiently decrease the planet's semi- major axis when 
either the so-called stochastic instability (see e.g. Kochanek (1992), Kosovichev 
k Novikov (1992), Mardling (1995)a,b, Ivanov & Papaloizou (2004), (2007)) 
operates in the dynamical system consisting of the orbit and tidally perturbed 
normal modes of the planet and the star, or when the normal modes are 
able to dissipate the energy transferred from the orbit during a time of order 
of or smaller than one orbital period. As we discuss below both processes 
can be effective when the semi-major axis of the orbit is sufficiently large, of 
the order of several au. Therefore, in general, the time t ev characterises only 
the orbital evolution at such scales. At smaller scales, nonlinear effects not 
considered here, such as wave breaking, or frictional processes acting on quasi- 
static tides can further decrease a thus completing the process of circularisation 
and leading to the formation of close-in giant planets with small P b s and 
eccentricities (e.g. Barker & Ogilvie 2009 and references therein). Provided 
that typical time scale of orbital circularisation due to these processes is smaller 
than t ev , this characteristic time scale may be considered as giving an overall 
time for the planet to circularise from an initially large semi-major axis a. In 
the opposite limit the overall time scale is obviously given by the time scale 
associated with the later evolution. 

In order to evaluate t ev we have to know the dependence of the stellar 
angular velocity, Q, on time t. In what follows we use the simplest possible 

1 As far as we are aware, Lai was the first to point out that energy transfer during a flyby 
on a retrograde orbit is enhanced compared to the prograde case. In this work the orbital 
plane was taken to coincide with the equatorial plane. Here we allow for the general case of 
arbitrary mutual orientations of these planes. 
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form of such a dependence assuming the so-called the Skumanich (1972) law 
for stellar rotation Q cx i -1 / 2 , where the coefficient of proportionally is chosen 
in such a way that the star has a period of rotation of one month at an age of 
5Gyr corresponding to that of the Sun. 

We find t ev to be different for prograde and retrograde orbits, with the 
difference increasing with planet mass. The ratio of the evolution time scale 
corresponding to a prograde orbit to that corresponding to a retrograde orbit 
is of order 1.5 — 2 for a one Jupiter mass planet with a typical P b s ~ 4d, and it 
is of order 5 for a five Jupiter mass planet with the same period. This difference 
gets more pronounced in a situation where tides exerted in the planet are not 
taken into account. Since conditions determining the long term effectiveness of 
the transfer of energy from the orbit to modes of pulsation can be different for 
the planet and star, such a situation may indeed occur. The results reported 
in this paper may be directly used in numerical experiments involving planet- 
planet scattering which also include tidal interactions between planets and the 
host star. 

The plan of the paper is as follows. We develop an analytic theory for cal- 
culating the tidal response of a slowly rotating star due to close the passage 
of a planetary perturber in a parabolic orbit in section [5J The orbital plane 
is allowed to have an arbitrary inclination to the stellar equatorial plane and 
corrections to stellar normal mode frequencies correct to first order in the stel- 
lar angular velocity are taken into account. We give expressions for the energy 
transferred to the star in section [3] and the analogous expressions for the trans- 
ferred angular momentum in Appendix B. In section 0] we describe numerical 
solutions of the linear tidal encounter problem. Following our previous work 
(Papaloizou & Ivanov 2010) we take account of Coriolis forces, but not cen- 
trifugal forces so that the tidally undisturbed model was spherical. However, 
apart from this, the stellar angular velocity may take any value. We perform 
calculations for a polytrope of index n = 3 and unlike in our previous work we 
include a treatment of the self-gravity of the star. 

Results for the energy and angular momentum transferred to the star as 
a consequence of prograde encounters, retrograde encounters and encounters 
with orbital plane perpendicular to the equatorial plane are then given in sec- 
tion [5] We then go on to compare numerical and analytic results in section [5] 
In both approaches the tidal transfers are significantly larger for retrograde 
encounters as compared to prograde encounters. We use our results to esti- 
mate the initial time scale for tidal circularisation for an orbit with an initial 
semi- major axis of 10au that circularises in a short period orbit in section I5T21 
for planets of one and five Jupiter masses. We give estimates with and with- 
out contributions arising from tides induced in the planet obtained from our 
previous work (eg. Ivanov & Papaloizou 2007). The central star is taken to be 
a solar mass star with angular velocity determined by the Skumanich (1972) 
law. We find that tides induced in the star and planet may separately account 
for the initial circularisation of orbits down to periods of a few days. Finally 
in section [7] we summarize and our discuss our results and their application to 
the tidal circularisation of prograde and retrograde Hot Jupiters. 
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2 The parabolic encounter of a planet with a rotating star: Analytic 
development for slow rotation and orbits of arbitrary inclination 

We develop the theory of the tidal interaction of a planet on a parabolic orbit 
with a central star. When the star is not rotating the tidal interaction excites 
p, f and g modes. A slow rotation approximation is then adopted in which the 
first order corrections to the mode frequencies due to rotation is taken into 
account in calculating the tidal response for an orbit with arbitrary inclination. 



2.1 Coordinate systems 

We calculate the energy transfer occurring as a result of the action of dynamic 
tides induced in a rotating star by a planet moving in a highly eccentric orbit 
in a plane that is inclined with respect to the stellar equatorial plane. The 
eccentricity is taken to be very close to unity so that to a good approximation 
the planet may be viewed as undergoing a sequence of parabolic encounters 
with the central star. 

For such a configuration, there are two natural Cartesian coordinate sys- 
tems with origin at the centre of the star. The first, (X' ,Y' , Z'), is such that 
the (X',Y') plane coincides with the orbital plane. The second, (X, Y, Z), is 
such that the (X, Y) plane coincides with the equatorial plane of the star. It 
is assumed that the the line of nodes which coincides with the Y axis makes 
an angle 7 with the Y' axis and lies in the interval [0,27r]. The direction of 
orbital motion is assumed to be correspond to a rotation about the Z' axis 
when viewed in the (X 1 , Y') plane with the apsidal line coinciding with the X 1 
axis and pericentre such that X 1 > 0. Thus the (X', Y') plane coincides with 
the orbital plane with the orbit being symmetric with respect to the X' axis. 
The (X, Y) plane coincides with the star's equatorial plane with the Z axis 
directed along the star's rotation axis. The line of intersection of the (X',Y') 
plane and the (X, Y) plane is chosen to define the Y axis. 

The angle j3 between the Z') and Z axes is called the inclination angle. Its 
values lie in the range [0, n]. Clearly, j3 = 7r corresponds to a situation when the 
planet's orbit lies in the equatorial plane of the star but the direction of orbital 
motion is opposite to the direction of stellar rotation. It is easy to see that the 
mutual orientation of the coordinate systems defined in this way corresponds 
to the angles j3 and 7 being the usual Euler angles that are adopted to define 
the relative orientation of two coordinate system^. Additionally, we make 
use of two spherical polar coordinate systems with spherical angles (#',<//) 
and (#,</>) respectively. The pairs of angles (#',</>') and (0,(f>) are defined in 



2 Note that the definition of the Euler angles correspond to the so-called zyz convention. 
Another way of defining the angles would be to perform the first rotation about the X axis 
and the second one about the Z axis. The so-defined Euler angles j3 and 7 correspond to 
the frequently used 212 convention. They are related to the angles used in this paper as 
P = j3 and 7 = 7 — tt/2. 
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the usual way, but with respect to the coordinate systems (X' ,Y' , Z') and 
(X, Y, Z), respectively. 



2.2 The tidal potential 

The quantity that determines the energy and angular momentum transferred 
between planet and star is the tidal potential, U. When the (X', Y', Z') system 
is used the tidal potential has the standard form (see e.g. Press & Tcukolsky 
1977, hereafter PT) 

r*/u>2 

U = W E W m e- im ^Y 2 , m (6', </>'), (1) 

where G is the gravitational constant, M is the planet's mass and r is the 
spherical polar radial coordinate. R(t) and <£(i) are the distance between the 
planet and the star and the angle between the position vector of the planet 
and the X axis measured in the orbital plane respectively. The time t = 
corresponds to a passage through periastron by the planet. The prime above 
the summation sign implies hereafter that the summation is performed over 

in = 0, i2 only. For \m\ = 2, W m = x figr and W = -J%, with Y Lm being 



10 

the usual spherical functions. 

When we consider the tidal potential in the (X', Y', Z') coordinate system 
with orientation determined by the stellar rotation axis, its form is modified. In 
order to find the new form of U we can use the transformation law of spherical 
functions under a rotation of the coordinate system. This may be written 

Y 2 , m (6', 0') = e-^ d$ n Y 2 , n (8, <f>), (2) 

n 

(2) 

where d n ;m are the so-called Wigner d functions, see e.g. Varshalovich et al 
(1989). Here the summation is performed over all integers, n, such that n £ 
[—2, 2]. Note that the d functions depend only on the inclination angle, /3, and 
satisfy the condition d> n _ m — {— \) n ~ m d_' n m . Their explicit forms are given 
in Appendix A. 

Substituting equation © in ([1]) we get 

C/ = r 2 ^A„Y 2 ,„(0,0), (3) 

n 

where 

Note that from the requirement that U is a real quantity and the properties 
of the spherical functions that F 2 *n — (~l) n ^2,-n, where (*) stands for the 
complex conjugate it follows that 

Al = (-l)M_ n . (5) 
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In order to evaluate the effects of a single encounter of the planet with the 
star, we introduce the Fourier transform in time of any quantity of interest, 
say F(t), F(a), by the relations 



/-r-oo i r-foo 

daexp(-iat)F(a), F(a) = — / dtexp(iat)F(t). (6) 

From equations ([5]) and (JSJ) it follows that the Fourier transforms of the am- 
plitudes A n obey the symmetry relations 

A:(-a) = (-l) n i-„(cr). (7) 

As follows from equations ^ and ([3]) the time dependence of the tidal 
potential is determined by the quantities I m = {R m i n / R(t)) 3 e~ lm ' p ^ t \ where 
Rmin is the periastron distance. Assuming that the orbit is formally parabolic, 
R{t) and <P(t) can be obtained from the parametric representations 

-^r— = l + * 2 (8) 

<P = 2t&n~ 1 x (9) 
n orb t = y/2 ( x + - oo < t < oo (10) 



Using these the Fourier transforms of these quantities, I m , can be readily 
obtained in the form 

Im = -^-h,-m{vn-\), (11) 
7TJ torb 

where 



CM 

n orb = J(l + q)^ (12) 



it a characteristic angular frequency of periastron passage, M* is the mass of 
the star, q = MjM* 1 is the mass ratio, and 

f°° dx 

h,m{y)= I 7z ttttt cos(V2y(x + x 3 /3) + 2marctana;), (13) 

Jo {l + x 2 ) 2 

see e.g. Pljf|. Using equations dU) and the definition of I m we obtain as explicit 
expression for A n (a), 

M°) = if Z M °n, a n = J2w m e- im idW m I 2 ^ m . (14) 



3 Note a misprint in Ivanov & Papaloizou (2007). In their equation (61) equivalent to (1 1 3 I t 
there should be I2, m on the left hand side instead of 72, -m as in the paper. 
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2.3 Equations of motion 

As has been discussed in our previous papers (Ivanov & Papaloizou 2004, 
Papaloizou&Ivanov 2005, Ivanov&Papaloizou 2007, Ivanov & Papaloizou 2010, 
Papaloizou&Ivanov 2010) for the case of coplanar encounters a perturber mov- 
ing on eccentric orbit excites normal modes of a rotating body after the pe- 
riastron passage. We must, accordingly, solve equations of motion of stellar 
perturbations associated with the modes under the forcing determined by the 
tidal potential in order to find the transfers of energy and angular momentum 
from the orbit to the modes. As we discuss above in this paper we would like 
to relax the assumption of coplanarity of the orbital and equatorial planes 
treating, however, the effect of stellar rotation in a simplest possible way and 
assuming that the angular frequency of the star, Q 1 is much smaller that an- 
gular frequencies of the normal modes, io nm , mainly excited in course of the 
tidal encounter, fi <C oj nm . In this case the normal modes frequencies are close 
to the ones of non-rotating star, uj a , where the index a labels different types 
of normal modes of a non-rotating star (e.g. p-modes, g-modes and the fun- 
damental mode). In the linear approximation in small parameter J7/w a only 
the modes frequencies acquire corrections of order of i? due to rotation while 
other quantities remain the same as in the non-rotating case. In the same 
approximation the total energy and angular momentum transfers can be rep- 
resented as a direct sum of partial contributions corresponding to modes of 
different type, and, therefore, we derive below these expressions for a mode 
of one particular type, omitting, for simplicity, the index a below, thus e.g. 
ui = uj a , etc.. 

The stellar perturbations can be described in terms of the Lagrangian 
displacement vector, £(i, r). Since the tidal potential is represented as a sum 
of the spherical harmonics Y/ jTl having I = 2 and —2 < n < 2 (see equation (J3])) 
in order to solve the response problem for a mode having in the non-rotating 
limit a frequency u> and I — 2 we should consider the vector £(£, r) as a linear 
combination of the eigen vectors with the same values of n, 

£ = £M*)Ur), ( 15 ) 

n 

where 

Ur) - Ur)Y 2 ,ne r + £ s (r)(rVY 2 ,„), (16) 

e r is unit vector in the radial direction, the functions £ r and t; s can be found 
by solving the eigenvalue problem for a non-rotating star, see e.g. PT, Lee &: 
Ostriker 1986. It is assumed hereafter that the eigen-vectors of the star are 
normalised by the condition (PT) 

d 3 x P (C ■ €») = [ ' r 2 drp{g + 6£ s 2 ) = 1, (17) 
Jo 

where integration is performed over the volume of the star, p is the stellar 
density. 
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Equations of motion for the coefficients b(t) can be found e.g. in Ivanov & 
Papaloizou (2004) as 

b n + w%b n + 2iu x n b n = S, with S = / d 3 xp(C n ■ VE7), (18) 



where for a normal mode with time dependence through a factor exp(— iwt), 
we set to = loq + w^, with ujq being the eigenfrequency for the non rotating 
star, and the frequency correction due to rotation of the star is w„, which can 
be represented as (see e.g. Christensen-Dalsgaard 1998) 

w* = nf3 r f2, (19) 

the values of the coefficient j3 r ~ 1 were calculated by Saio (1981) for a number 
of I = 2 normal modes of a n — 3 polytropic star. 

Substituting equation (O in the expression for the source term S in (fT5|) 
we get 

f R, 

S = A n Q, Q= / r 3 drp(2^ + 6^ s ). (20) 
Jo 

The overlap integrals Q have been calculated by PT and Lee & Ostriker 1986 
for a number of polytropic stellar models. 

From equation (JT5J) it follows that the Fourier transform of the amplitude 
6„, b ni can be expressed through the Fourier transform of the amplitude A n 
defined in equation A n , as 

k = ^r\, D(a)=u%-(a + w-wl) 2 , (21) 

where we add a small imaginary correction tocr, a — > a + iv, v > according 
to the Landau prescription and neglected (w^) 2 as a small quantity. 



3 The energy transfer 

During the planet's flyby of the star the energy is deposited in the star at a 
rate determined by work done by the tidal force 

E = J d 3 xp(£-VU), (22) 

and, accordingly, the total energy transfer, AE, is obtained by integration of 
(|2"21 over time 

dt / d 3 ccp(i • V£T), (23) 

-oo J 

where the limits of integration are formally extended to ±oo since the integral 
over time in (I23[) is highly peaked near t = corresponding to the periastron 
passage. 
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We substitute equations (|T5j) and (HI))) in ([2"5")) and integrate over the 
volume of the star using the known orthogonality properties of the spherical 
functions and relations ([5]) to get 

/+oo 
dtb n A* n , (24) 
-°° 

where we recall that the overlap integral Q is given by (|20|) . Now we express 
b n and A n in (|2"B1 in terms of their respective Fourier transforms using © 
and integrate the result over time with help of the well known property that 
/ dte-^-" ') = 2ir5{cr - a') to obtain 



AE 



-2™Q 2 ]T I™ dJ-^f- (25) 



where -D(cr) is defined in equation (|2"Tj) . 

A simple analysis shows that in the limit of v — > the integral in (1231) is 
mainly determined by contributions close to the poles of the expression under 
the integral, where D(a) « and, accordingly, 

tr«w±=±w +w* (26) 

Taking into account only these contributions in (|25l) we get 

Z\£ = (ttQ) 2 ^ (l + ^ |i„(c+)| 2 + (l " ^) IA.«)| 2 . (27) 

It is convenient to use only positive values of the resonant frequencies in (|27[) . 
Taking into account that uo 1 _ n — —co^ and, accordingly uj~ — — w_„, and using 
equation (0 we have 

= („0) 3 £ (l + ^) |i„Ht)| 2 + (l- ^) |I-.,( U ;„)| 2 . (28) 



Now we substitute equation (|14l) in (|28p to obtain 
where 

e = |a H| 2 + f 1 + ^ |«iK)| 2 + f 1 - ^ ) lo-ifwij) 



1 + _£ I loaK)!" + I 1 - ^ I |a_ a («± 2 )r, (30) 

and we take into account that the terms with opposite signs of n give the same 
contribution in the summation series in (I28[) . Using equation (|14l) we can also 
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obtain explicit expressions for the quantities |a n (w)| 2 entering in (|30j) in the 
form 



|a n (w)| 2 = Wl{d\ 2 l\_ 2 + d 2 n _ 2 ll 2 + 2d n ^d n -2h,-2h,2 cos 4rf) 
+ 2W 2 W d nfi I 2fi cos2 7 (d„ !2 / 2 ,_ 2 + d„,_ 2 / 2 ,2) + W$d 2 nfi ll fi , (31) 



where all Ii^ n are assumed to be functions of u> and we set, for simplic- 

(2) 

ity, d nt m = dn'm- One can easily check that when /3 = and, accordingly, 
dn.m = $n.m the expression ([29]) is equivalent to the expression (59) of Ivanov 
& Papaloizou (2004). 

It is instructive to express natural units in (|29[) . Taking into account that 
according to our normalisation condition (|17[) the overlap integral Q can be 

1/2 ~ 

expressed as Q = i?»M* Q, where i?* is the radius of the star and Q is dimen- 
sionless. Let us also introduce natural energy and frequency scales associated 
with the star, E* — GM 2 /i?* and 12* = y/ GM*/Rl, respectively, and use the 
Press Teukolsky parameter rj characterising a strength of the tidal interaction, 



With help of these new variables the expression (|2"9")1 can be rewritten in the 
form 



4 Numerical solutions of linear tidal problem for orbital planes 
inclined to the equatorial plane of the star 

Here we describe numerical solutions for tidal encounters of a planet with a star 
when the orbit is not coplanar with its equatorial plane. We obtain the energy 
and Z component of the angular momentum exchanged due to tides acting 
on the star. This in turn can be related to the rate of orbital circularisation 
starting from high eccentricity (Papaloizou & Ivanov 2007). 

The numerical procedures follow closely those described in Papaloizou & 
Ivanov (2010) that were applied to tidal encounters with a polytrope of index 
n = 1 but with appropriate modifications to consider a polytrope with n — 3. 
These are described below. Only coreless models are considered here. As in 
Papaloizou & Ivanov (2010) grid resolutions of 200 x 200 and 400 x 400 have 
been considered. 




(32) 
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(a) (b) 

Fig. 1 a) The energy, in units of E-* , transferred to the non rotating star as a result of a 
parabolic encounter with r\ = 4\/2 (upper solid curve) and rj = 8 (lower dotted curve). The 
mass ratio q = 10 — 3 . 

b) The z component of the angular momentum, in units of J* = M* \J GM-* R* , transferred 
to the non rotating star for r\ = 4\/2 (upper solid curve) and r) = 8 (lower dotted curve). 
The mass ratio q = 10~ 3 . 



4.1 A simple stellar model 

We adopt a simple model corresponding to a spherically symmetric polytrope 
of index n = 3. It is in hydrostatic equilibrium such that 

i£ = * (34) 
p dr dr 

where p is the density, P = Kp 1+1 / n is the pressure, with K being the poly- 
tropic constant and ip is the- gravitational potential arising from the stellar 
matter. This satisfies the Poisson equation 



(35) 
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4.2 Linearised equations of motion 

We assume that the star rotates with uniform angular velocity fl directed 
parallel to the Z axis. The hydrodynamic equations for adiabatic perturbations 
induced by an orbiting planet in a frame corotating with the star are 

Qv p 1 /i f 

— + 2S2 x v = V(P'/P 1/7 ) - &.c4 v e r + VU + V^' + — , (36) 

at p p 

where the Brunt Vaisala frequency ujbv is given by 

a _ ldP fldp 1 dP\ 

v = (w r ,ug,W0) is the Eulerian velocity perturbation, P' is the Eulerian pres- 
sure perturbation, f„ is the viscous or diffusive force per unit volume, ip 1 is the 
perturbation of the gravitational potential due to the star and as indicated 
above U is the external forcing tidal potential. The adiabatic exponent 7 is 
taken to be constant and equal to 5/3. The velocity perturbation v is related 
to £ = {£ r ,£e ,£</>), the Lagrangian displacement vector, through 

The linearised continuity equation gives 

P = -V • (p€), (39) 

where p' is the density perturbation. This together with the adiabatic condition 
gives the Eulerian pressure perturbation as 

dP 

P' = -£ r — - 7 PV-(pO- (40) 

The gravitational potential perturbation tp' satisfies the linearised Poisson 
equation 

VV = -4ttGp'. (41) 

Equation Q41]i may be solved by expanding the right hand side in a series 
of spherical functions. We simplify the problem by noting that we consider 
solutions the tidal problem by combining solutions of equation (|36l) where the 
angular dependence of U is through a single spherical function. To evaluate 
i/j' we retain only the spherical function associated with the forcing in the 
decomposition of p' . This is an exact procedure in the non rotating case and 
should be a reasonable approximation in the rotating case as the potential 
perturbations associated with high order spherical functions are expected to 
be small. It amounts to retaining only the quadrupole components of the 
gravitational potential for the problem on hand, a procedure that has been 
used successfully in problems of stellar dynamics (eg. Allen et al. 1990). 
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Note also hat the centrifugal term is absent in equation (1361) being formally 
incorporated into the potential governing the static equilibrium of the unper- 
turbed star and there is no unperturbed motion in the rotating frame. As in 
our previous work (eg. Papaloizou & Pringle 1981, Papaloizou & Ivanov 2010) 
we shall neglect centrifugal distortion of the basic equilibrium which enables 
us to adopt a spherically symmetric unperturbed model and density distribu- 
tion. We consider prograde encounters, retrograde encounters and encounters 
for which the orbital plane is at right angles to the stellar equatorial plane. 



4.3 Tidal forcing potential for prograde and retrograde encounters 



For prograde encounters the general expression ([3]) for the tidal potential re 
duces to the specific form 

GM r 2 



U = K 



3GM / r \ 2 
4i? \R 



sin 2 (6») exp(2i((t) - $) 



AR R 2 



(3cos 2 0- 1) 



(42) 



where 1Z indicates that the real part is to be taken and R(t) and <P(t) are given 
by © - (fTT)|) . The same expression may be used for a retrograde encounter by 
making the replacement <P — > —<P. 



4.4 Tidal forcing potential for encounters with j3 = it/ 2 

In this case when pericentre is in the equatorial plane we may write 

U = — — {—) sin 2 9 cos 2 <2> exp(2#) (43) 

4it \RJ 

-6i sin 9 cos 9 sin 2<£exp(ic/>) + (3cos 2 6>- l)(3sin 2 ^- 1)] 
and when it is on the rotation axis we may write 

GM ( r \ 2 r 99 , , . 

U = — j K [-3 sin 2 9 sin 2 <P exp(2i0) (44) 

+6i sin 9 cos 9 sin 2$ exp(i^) + (3 cos 2 9 - 1)(3 cos 2 $ - 1)] . 

We remark that changing the direction of motion along the orbit by the re- 
placement <P — > — <P does not change the energy transfer or the change in 
the Z component of angular momentum. This corresponds to encounters with 
(3 = ir/2 and (3 — 3tt/2 giving the same results for these quantities. The above 
tidal potentials consist of a sum of terms with different azimuthal mode num- 
ber m. These can be considered independently and the energy and angular 
momentum transferred as a result of the action of each term evaluated as the 
canonical energy and angular momentum obtained after the encounter is over 
(see below). The results can then be added to get the total. As in Papaloizou 
& Ivanov (2010) the encounters are started with the perturber at eight times 
the pericentre distance with R and <P determined from © - (fl"U|) . It can be 
seen that the results are unaffected by applying an arbitrary complex phase 
shift to each individual forcing term. 
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Fig. 2 a) The energy, in units of £*, transferred to the star for rj = 4\/2 and q = 
10 -3 plotted as a function of fi/fi*. Negative values of the latter correspond to retrograde 
rotation. The upper dashed curve spanning (—0.25, 0.25) shows the energy transferred when 
the orbital plane is at right angles to the stellar equatorial plane where pericentre is located. 
The lower dot dashed curve spanning (—0.25, 0.25) is for the same parameters except that 
pericentre is located on the stellar rotation axis, b) The angular momentum, in units of 
J* = M* V GM t Rt, , transferred to the star. The initially uppermost dashed curve spanning 
(—0.25, 0.25) shows the z component of the angular momentum transferred when the orbital 
plane is at right angles to the stellar equatorial plane where pericentre being is located. 
The initially lowermost dot dashed curve spanning (—0.25, 0.25) corresponds to the same 
parameters except that pericentre is located on the stellar rotation axis, c) As in a) but for 
r] = 8. d) As in b) but for rj = 8. 
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4.5 Addition of diffusive effects 

In order to avoid potential numerical problems arising from the excitation 
the a dense or continuous spectrum of normal modes that may be associated 
with a rotating star, we follow Papaloizou & Ivanov (2010) and incorporate a 
diffusive force per unit mass, f„, in the simulations. For models with non zero 
rotation this is taken to be of the form 

{ - — (pv—\ + — — ( pv& m6—\ (45) 
v dr \ dr J r 2 d/j, \ dp, J 



where the diffusion coefficient or effective kinematic viscosity v — vq\/GM*R* 
was taken to be a constant. For simulations with a grid resolution of 200 x 200 
we adopted v$ = 1.3 x 10 -6 and for a grid resolution of 400 x 400 we adopted 
vq = 6.5 x 10~ 7 . For non rotating models vq = 0. We remark that energy and 
angular momentum transfers measured by evaluating the canonical energy 
and angular momentum just after the encounter has completed are robust to 
changes in numerical resolution and diffusivity. 

4-5.1 The canonical energy 

When the Lagrangian pressure perturbation vanishes at the outer boundary, 
the canonical energy (defined in the rotating frame) appropriate to forcing 
through a Fourier mode with azimuthal mode number, to, may be written as 



E c = 0.25(1 + S mfi ) 



P (\v\ 2 +UJ% v \^\ 2 + \PT/( 7 P))dT 



/^G~ dT ~ J A \Zr\ 2 (dP/alr)dS 



(46) 



where <5 TOi o) is the Kronecker S. The canonical energy contains contributions 
from the pressure contribution, buoyancy and self-gravity. The last two were 
not included in Papaloizou & Ivanov (2010) where a polytrope with n = 1 was 
considered using the Cowling approximation. The first volume integral is taken 
over the volume of the star V, while the second giving the gravitational energy 
is taken over all space. The surface integral is taken over the surface area 
A. The contribution of the surface term formally vanishes when the density 
at the surface is zero. When the density at the surface is relatively small as 
in our model, this term gives a negligible contribution. When tidal forcing 
operates, the time rate of E c gives the rate of energy transfer to the star. 
After the encounter tidal forcing tails away and E c is conserved in the absence 
of dissipation. In the presence of dissipation it decays with time. 

Let us stress the difference between the energy defined in the rotating frame 
and in the inertial frame. The energy transfer in the rotating frame AE r = E c 
is related to the energy transfer in the inertial frame (see equation (j33|) and 
the next Section ) as 

AE = AE r + QAJ Z , (47) 
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where J z is the transfer of zth component of angular momentum, which may be 
equated to the canonical angular momentum, J c , contained in the perturbed 
modes, see below. While it is more convenient to operate with AE r in our 
numerical work, the analytical work as well as different astrophysical applica- 
tions are concerned directly with AE. Therefore, we illustrate the behaviour 
of AE r in this Section when discussing our numerical results and AE in the 
next Section where numerical and analytic results are compared. 



4-5.2 The canonical Angular momentum 



The corresponding canonical angular momentum is given by 



0.5m / p (C r v r + $v + Qv* + 2rt(&.Q sin 9 + cos 9)) dr 



■I, = -2m 

(48) 

where Tm indicates that the imaginary part is to be taken. Here we recall that 
the state variables are complex. This behaves in a manner analogous to the 
canonical energy, but as applied to the total angular momentum content of 
the star. 



5 Numerical results 



5.1 The energy and angular momentum transfer for a close encounter 



The energy, in units of E*, transferred to a non rotating star during parabolic 
encounters with a planet with mass ratio q = 1CP 3 for r\ = 4\/2 and n = 8 
is shown as a function of time in Fig. 1(a) These results are for the m = 2 
tide which is the dominant one (Papaloizou & Ivanov 2004, 2010). The corre- 
sponding angular momentum transferred to the star is plotted as a function of 
time in Fig. 1(b) The energy and angular momentum transferred level of at 
constant values after the tidal interaction is essentially complete, the charac- 
teristic time associated with the encounter being about 1 day. This corresponds 
to the excitation of the normal modes associated with the rotating star which 
in this case correspond to /, p and g modes. It is found that g modes up to 
order 9 may be significant (see Lee & Ostriker 1986). In practice we measure 
transferred quantities at t = 0.75d when n = 4\/2 and at t = lAd when n = 8. 
Our results for the final energies transferred to the star agree with theirs to 
about the 10% level. 



5.2 Prograde and retrograde encounters with a rotating star 



The energy, in units of E* 
parabolic encounter with n 



transferred to the star as on completion of a 
= 4y2 and q — 10~ 3 is plotted as a function of 



4?/J2* in Fig. 2(a) Negative values of QjQ* correspond to retrograde rotation. 
The angular momentum, in units of J* = M*^JGM*R*, that is transferred is 
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Fig. 3 a) The energy, in units of E*, transferred to the star for rj = &\/2,q = 10~ 3 and 
J7/J7* = 0.24 (uppermost solid curve). The orbital plane was at right angles to the stellar 
equatorial plane with pericentre located in the latter plane. The initially lowermost dashed 
curve shows the contribution from terms in the forcing potential with m = 2. Similarly the 
initially next lowermost dotted curve corresponds to m = 1 and the ultimately lowermost 
dot-dashed curve to m = 0. b) The z component of the angular momentum , in units of J* , 
transferred to the star (lowermost solid curve). The ultimately intermediate dashed curve 
shows the contribution from terms in the forcing potential with m = 2. The ultimately 
uppermost dotted curve corresponds to m = 1. c) As in b) but for rj = 8 and Q/Q-+ = 0.12. 
d) As in c) but for r) = 8 and i7/J7* = 0.12. 
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plotted as a function of fl/fl„ in Fig. 2(b) In Fig. 2(c) and 2(d) the corre- 
sponding plots are given for rj = 8. The energy transferred decreases as fl is 
increased from zero until a minimum value is attained, after which it increases 
with fl again. On the other hand the energy transferred increases with in- 
creasing retrograde rotation attaining values up to more than three orders of 
magnitude larger than the minimum one for fl/ fl* ~ —0.6. when rj = 8. Thus 
circularisation from high eccentricity as a result of stellar tides can be expected 
to be much more effective for significant retrograde rotation. For r/ = 4v2 the 
minimum energy transfer occurs for fl/fl# ~ 0.4 while for rj = 8 this occurs 
for fl/fl* ~ 0.25. The actual angular velocity at pericentre flp/fl„ = v2/r?, 
so that fl p /fl* = 0.25 and 0.18 for these two cases respectively. Accordingly 
this minimum transfer occurs when fl ~ 1.5fl p , which is approximately the 
angular velocity for which the angular momentum transfer reverses sign, this 
transition occurring when fl ~ fl p , (see Figs. 2(b) and 2(d) I. This character- 
istic behaviour was also noted in the case of a polytrope with index n = 1 and 
a barotropic equation of state ( Ivanov & Papaloizou 2004, 2007). 

We have also considered encounters for which the orbital plane and the 
stellar equatorial plane are perpendicular. Cases for which pericentre was in 
the equatorial plane and on the rotation axis have been considered. The energy 



transferred to the star for some of these cases is plotted in Figs. 2(a) and 2(c^ 



The component of angular momentum in the z direction transferred to the star 
is plotted in Figs. 2(b) and |2(d)| The energy transferred in these cases is also 
seen to significantly exceed the minimum that is transferred in the prograde 
case. 

To look at such encounters more closely, the energy, in units of E*, trans- 
ferred to the star as a result of a parabolic encounter with 77 = 4y2 and 
q = 10~ 3 for fl/fl* = 0.24 when the orbital plane was at right angles to the 
stellar equatorial plane where pericentre was located is plotted as a function 
of time in Fig. 3(a) The z component of angular momentum transferred is 
plotted in Fig. 3(b) The contributions from the terms in the forcing potential 
with m = 2, m = 1 and m = are indicated. Clearly the last of these does not 
affect the z component of angular momentum. Note that significant contribu- 
tions come from m = 1 as well as m = 2. The corresponding results for 77 = 8 
and fl/fl* = 0.12 are shown in Figs. 3(c) and|3(d)] 



6 Comparison of numerical and analytic results and circularisation 
time scales 

6.1 Numerical versus analytic results 

In order to calculate the analytical value of the energy transfer, AE, we should 
know the values of the mode eigen frequencies, the overlap integrals and the 
frequency splitting coefficients /3 r . We take into account the /-mode, 5 low 
order p-modes and 14 g-modes, which give the most important contribution 
to the energy transfer. The eigen frequencies and the overlap integral are taken 
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a) 




b) 



Fig. 4 a) Analytic and numerical dependencies of the energy transferred to the star in units 
of the characteristic energy E* (E* = 3.8- 10 48 ergs for a star of solar mass and radius) on the 
angular velocity of the star in units of the characteristic frequency i?* = 6.23 ■ 10 _4 s _1 
for a star of solar mass and radius), for r) = 4\/2- See the text for a description of the 
different curves and symbols, b) Same as a) but for i) = 8. 
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from the papers Lee & Ostriker (1986) and Ray et al (1987). The coefficients 
(3 r for the /, p and several low order g-modes are taken from the paper Saio 
(1981). For the rest of g-modes we use the WKBJ results that j3 r ~ 1 — X(T+T) > 
where L = 2 in our case, see e.g. Christensen-Dalsgaard (1998). 

In Figs. 4(a) and |4(b)j we show the dependencies of the energy transferred 
to the star as a result of the planet flyby on the angular velocity of the star, for 
■q = 4^2 and rj = 8, respectiveljQ. The solid, dashed and dot dashed curves 
represent the analytic results while the symbols indicate the corresponding 
numerical results. The solid curve and circles correspond to the situation when 
the equatorial and orbital planes coincide (/3 = 0). The dashed curve and 
squares correspond to an orbit with normal to the orbital plane perpendicular 
to the rotational axis and the direction to pericentre laying in the equatorial 
plane (/3 = n/2 and 7 = t/2). The dot dashed curve and diamonds are 
calculated for the same f3 = ir/2 but for the direction to pericentre being 
aligned with the rotational axis, and, accordingly, 7 = n. Note the one can 
either consider negative values of ft as in Figs. 3(a)|3(d)| or a change of the 
inclination angle such that /3 — > /3 + it keeping the value of ft positive. Both 
approaches are physically equivalent and give the same results. In what follows 
we describe the results corresponding to positive (negative ) values of ft as 
prograde (retrograde ) and the results corresponding to/? = 0( j 8 = 7r/2) 
as planar (perpendicular). As seen from Figs. 4(a) and |4(b)"| in the planar 
case the analytic and numerical results corresponding to positive rotational 
frequencies are in good agreement, while the results corresponding to negative 
frequencies deviate, leading to systematically larger analytic values of AE as 
compared to the numerical determined ones. The deviation is of the order or 
smaller than 50% for — ft) ft* < 0.1. For the largest value of — ft/ ft* w 0.25 
shown, it is about a factor of two for 7/ = 4-^2 and about a factor of 2.5 for 
rj = 4\/2. In the perpendicular cases both analytic and numerical results show 
much weaker dependence on fl. The analytic curves again correspond to larger- 
values of AE, with results deviating most significantly from the numerical ones 
at large values of The degree of deviation behaves in a similar manner 
to the planar case. Note that the analytic curves always grow with increasing 
\fl\ while the numerical curve corresponding to 7 = 7r/2 (7 = 0) is weakly 
increasing (decreasing ) with This weak growth (decrease) is probably 
within numerical errors and should be treated with caution. 

The origin of the deviation at large values of — ft/ fl* is unclear at present. 
It cannot be fully accounted for by the increasing importance of non-linear 
corrections to the rotational shift of eigen frequencies given by equation (fT9")) 
with increasing \Q\. We found that phcnomcnological non-linear corrections 
provided by Lai (1997) do not significantly affect our analytic results for the 
range of fl considered. However, for the values of r\ considered, the most sig- 
nificant contributions to the energy transfer are provided by g-modes of order 
of 9 — 10. The eigenfrequencies of these modes being ~ 0.5J?* are a factor of 



4 Let us recall that here we show the energy defined in the inertial frame, which is related 
to the energy defined in the rotating frame discussed in Section [5] through equation 1471 . 
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~ 2.5 larger than the maximum angular velocity considered. Thus they are 
approaching the so called inertial regime for which eigenfrequencies are less 
than 1 2 J? | in magnitude. It is possible that the overlap integrals become sig- 
nificantly modified in such a situation. This issue will be explored in future 
work. 



6.2 The time scale for the initial stage of tidal circularisation 
6.2.1 A basic approach 

In principle, the results reported above can be used in numerical experiments 
in which a system of gravitationally interacting planets is considered, it being 
assumed that once a planet has obtained a sufficiently small value of its or- 
bital angular momentum as a result of planet-planet gravitational scattering, 
it can lose orbital energy due to tidal interaction with the central star during 
successive periastron passages. If the effect of tidal interaction on changing 
the orbital energy becomes much larger than residual effects associated with 
planet-planet scattering, the planet is said to be tidally captured by the star. 
After being tidally captured, the planet interacts mainly with the central star 
due to dynamic and quasi-static tides exerted in itself and in the star. Under 
certain conditions discussed in the next section the orbital energy decreases, 
either systematically or on average, while the orbital period and orbital eccen- 
tricity also decrease. In this way it is possible to account for at least some of 
the so-called "Hot Jupiters" - giant planets with orbital periods ranging from 
a fraction of a day to several days. 

The central star can rotate, with its rotation axis inclined with respect to 
the normal to the orbital plane by an angle (3, after the tidal capture. This 
angle is determined by the particular sequence of planet-planet interactions 
that led to the tidal capture. Since retrograde tidal interactions that transfer 
energy to the star are stronger than prograde ones, while the angle j3 may be 
random or have a weak correlation with the direction of stellar rotation, this 
property of the tidal interaction may influence the observed relative number 
of Hot Jupiters on retrograde and prograde orbits. Thus it is important to 
estimate the potential magnitude of this tidal effect. 

As Hot Jupiters may have been formed when the star was young, in making 
such an estimate it is very important to account for the fact that the angular 
velocity of the star, Q, is a function of stellar age. This function can be a rather 
complicated and is expected to be different for stars with different masses and 
formation histories (e.g. Scholz 2009). In this paper we make only preliminary 
qualitative estimates of the importance of tides exerted in a rotating star on 
the formation of Hot Jupiters. Accordingly, we use the simplest plausible form 
of this function known as the Skumanich law (Skumanich 1972), which states 
that the angular velocity is inversely proportional to square root of the stellar 
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age. Thus we adopt the following expression for the stellar angular velocity: 



Here fl is normalised so that the star has a rotation period of one month at 
an age of t = 5 ■ 10 9 yr. 

In this paper we do not consider the planet-planet interaction stage, leaving 
this for a future work. Here we assume that tidal capture has happened and 
the planet interacts only with the star. It supposed to be initially on a highly 
eccentric ( essentially parabolic jl orbit with well defined values of the angles 
/3 and 7 and have a value of the orbital angular momentum corresponding 
to a given value of the parameter 77 defined in equation (|32p . The orbital 
energy, eccentricity and period are gradually changing their values due to 
tidal interactions on some characteristic evolution time scale, t ev , specified 
below. As discussed in e.g. Ivanov & Papaloizou (2004,2007), orbital evolution 
as a result of tides acting in the planet approximately conserve orbital angular 
momentum because the planet spins up (or slows down) by tides to a state 
of rotation corresponding to the so-called pseudo-synchronisation, where the 
angular momentum exchange between the orbit and pulsational modes of the 
planet is absent. This happens on a quite short time scale due to a small 
moment of inertia of the planet compared to that of the orbit. 

In the case of tides exerted in the star we have an opposite situation where 
the stellar moment of inertia is typically much larger than the orbital one. 
Therefore, the star cannot be brought by action of tides to the state of pseudo- 
synchronisation during the evolution of the system, and some amount of an- 
gular momentum is transferred from the orbit to the star in the course of it 
thus decreasing the value of the orbital angular momentum. We would like 
to estimate this value assuming that tides are effective all the way down to 
smallest values of the orbital eccentricity. Let us stress again that the dynamic 
tides discussed in this paper are effective only at sufficiently large values of 
the eccentricity e ~ 1 ( see the next section). However, it seems plausible to 
suppose that processes not considered here, such as the nonlinear breaking of 
excited waves, critical latitude phenomena in the planet, or friction, resulting 
from for example turbulence, acting on quasi-static tides could be important 
when the eccentricity is moderate and small, and that they can eventually 
bring the system to the state having a low eccentricity e ~ typical for the 
Hot Jupiters (e.g. Barker & Ogilvie 2009). 

It is well known that the mode energy, E m , is related to the mode angu- 
lar momentum, J m , as E m = uj m J m /m, see e.g. Friedman & Schutz (1977). 
Accordingly, the total amount of energy and angular momentum, E t and Jt, 

5 A planet with orbital parameters typical for the stage of tidal interaction right after 
the tidal capture would have its periastron distance of order of several (say, 4) stellar radii, 
while its semimajor distance would be of order of a few (say, 10) astronomical units (e.g. 
Nagasawa et al 2008). For such parameters orbital eccentricity e is of order of 0.998 and 
approximation of tidal encounters as happening on parabolic flybys is well justified. 




(49) 
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transferred from the orbit to the star obeys a similar relation, where, in gen- 
eral, we should sum over all modes of the star. In our case of a slowly rotating 
star only m — 2 o-modes of sufficiently high order give a significant contribu- 
tion to the sum. For our problem the corresponding eigenfrequencies of these 
modes can be roughly estimated as u> n ~ 0.5.!?*, and we obtain E t ~ f2* Jt/4. 
After the tidal circularisation has happened assuming that the tides induced 
in the star dominate over those in the planet we can equate the total mode 
energy E t to the orbital binding energy of the planet after the circularisation, 
Ej ~ GMM^/ (2a/), where a/ is the final value of the orbital semi-major axis, 
and we suppose that e ~ 0, and, accordingly, the orbital angular momentum 
Jf = MyJGM^af. From these expressions we obtain 

Jt -/i^ 3/2 



j f ~ 2 V a / J \/v ^ 

where in the last equality we use equation (|52")l setting q = there, and 
assuming that Jt/ Jf is small, and, therefore, the orbital angular momentum is 
approximately conserved during the tidal circularisation. In this case we can 

Set Rrnin 

= a//2 in ([211). As follows from equation (JCTJ) below, J t /J f < 0.1 for 



the Hot Jupiters with the orbital periods > 2.4 days. As seen from Figs. 5(a) 
and |5(bj] the stellar tides can dominate for even larger periods, accordingly the 
assumption of approximate conservation of the orbital angular momentum can 
then be justified.]. 

This fact also allows us to link a value of rj to a value of orbital period after 
the tidal circularisation process has been completed and the planet's orbit has 

became almost circular, P f, s = / (GM»), as 



'"'UsJ W p " (51) 

where P\ = P b s /^day. This equation relates an observable quantity, P h s , to 
the quantity determining the strength of the tidal interactions. 

For simplicity, it is assumed below in this section that we can simply add 
the values of energy transfer during successive periastron passages either due 
to sufficiently efficient dissipation of energy stored in the modes between pe- 
riastron passages (e.g. Kumar & Goodman 1996) or due to the operation 
of stochastic instability in the system, e.g. Kochanek (1992), Kosovichev & 
Novikov (1992), Mardling (1995)a,b, Ivanov & Papaloizou (2004), (2007). Note 
that in the latter case the characteristic time scales for orbital circularisation 
must be understood in some average sense. Some estimates relevant to a more 
realistic situation are discussed in the next section. 

Under this assumption Ivanov & Papaloizou (2004), (2007) derived an 
equation determining the rate of change of the orbital semimajor axis, a, in 
the form 



When the planetary tides dominate, equation (1 50 II clearly overestimates Jt/Jf. 
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aw 

where aio = a/(10au) and 




(52) 



7-10* — Vr, (53) 



where Mj, i?j are the Jupiter mass and radius, respectively, R p i is the planet's 
radius, and 

e DT = (AE + AE pl )(^^, (54) 

where AE is given by equation and Z\£" p / is the energy transfer to the 
modes exited in the planets . This is calculated in Ivanov & Papaloizou (2004), ( 
2007) for the value of rj given by (|5Tj) corresponding to a specified final orbital 
period. Note that both AE and AE p \ depend on time. In the case of AE this is 
through the dependence of the stellar angular velocity on time, assumed to be 
given by equation (|4T)]) . In the case of AE p i this is mainly through the evolution 
of the planet's radius on time. Following Ivanov & Papaloizou ( 2004), (2007) 
we assume in this paper that this evolution is the same as the evolution of 
an isolated planet in which the planet's radius shrinks with time. Thus, we 
neglect the possible influence of tidal heating on the planet's structure. 
Equation (|54"|) can be integrated to give 

( 1 f* dt ' \ 2 

a w (t) = a in 1 - -j= \ - , (55) 



2^i„ J Un tw(t') 

where ai n (ti n ) is the 'initial' value of the planet's semimajor axis (in units of 
10au) just after tidal capture has occurred. Let us define the evolution time 
scale as the time corresponding to aio being formally equal to zero. From (|55p 
we get 

dt' r- 
I 77a = 2 V« m , (56) 
tiolj ) 

which provides an implicit equation for determination of t ev . However, let us 
recall that the evolution time scale so defined describes the evolution of the 
semi-major axis at large scales, of the order of 1 — lOau or larger, under the 
assumption that the dynamical tides can effectively transfer the orbital energy 
to energy associated with normal modes and thermal energy of the star and 
the planet. At smaller scales the dynamic tides are clearly ineffective and other 
processes (e.g. due to so called quasi-static tides ) must be considered to explain 
final stages of circularisation of exoplanets, see the next section. 



7 Note that in Ivanov & Papaloizou (2004). (2007) there is a misprint in their equations 
defining tig. Namely, the value of t\o is proportional to there instead of the correct 

square root dependence of this ratio as in (153 D - However, in these papers only the case 
M„ = Mq was considered and this misprint has no influence on the results of the papers. 
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The results of the solution of this equation for aj„(ij n ) = 1 corresponding 
to an initial orbit with semi-major axis of lOau and a central solar mass are 
shown in Figs 5(a)| 5(b) for M = IMj and M = 5Mj, respectively. In both 



cases it is assumed that the planet is in a state of rotation corresponding to 
so-called pseudosynchronization, see Ivanov & Papaloizou (2004), (2007). We 
consider only the case 

of the planet's orbit being coplanar with the stellar equatorial plane but 
with both retrograde and prograde directions of orbital motion. As seen from 



Fig. 5(a) when tides are exerted both in the star and in the planet, there is a dif- 
ference in the circularisation time scale corresponding to the case of retrograde 
motion as compared to the cases of either a non-rotating star or of prograde 
motion. Since criteria for the dissipation of the energy in normal modes as well 
as conditions determining onset of the stochastic instability may be different 
for the planet and for the star and it may happen that only stellar tides are 
important for the circularisation problem under certain circumstances, we also 



plot the curves calculated assuming that AE p j — 0. As seen from Fig. 5(a) the 
stellar tides themselves can provide a tidal circularisation time scale smaller 
than a typical life time of a planetary system of ~ 5 • 10 9 yr for P or \, < 4 



days. In Fig. 5(c) we show the ratios of the evolution time scale corresponding 
to the prograde and retrograde cases to the case of a non-rotating star, for 
a planet with M — IMj. One can see from this figure that assuming that 
tides in the star and in the planet are both effective, there is practically no 
difference between the prograde and non-rotating cases, while the retrograde 
case has an evolution time scale 1.2 — 2 smaller, for final orbital periods in 
the range 3.5 to 4.5 days. When only the stellar tides are taken into account, 
the difference between the prograde, retrograde and non-rotating cases can 
be quite significant for smaller final orbital periods. Since AE grows with the 
mass of the planet the contribution of stellar tides is more important for the 



M = 5Mj case shown in Fig. 5(b) In this case t ev is practically determined 
by these tides for P or t larger than 2-2.5 days and the difference between the 
retrograde, prograde and non-rotating case is quite prominent. 

As mentioned above there is a difference between the numerically and an- 
alytically obtained values of AE, which is more significant for a sufficiently 
rapidly rotating star and retrograde orbits. Let us estimate the possible differ- 
ence in the evolution time scale associated with this assuming that t ev > 10 7 yr. 



From equation |g5J) it follows that Q{t > 10 7 yr) < 0.08/2*, and from Figs. [5(a) 



and 5(b) it follows that when J?//2* < 0.08 the ratio of analytically obtained 
values of AE to the numerical ones is of the order of or smaller than 1.5. In 
order to see the effect of tides exerted in the star being reduced by a factor 



1.5, we plot the additional dot dot dashed curve in Fig. 5(a) corresponding to 
the retrograde case with the energy transfer given by equation (|33|) reduced 
by a factor of 1.5. It is seen that the curve is very close to the short dashed 
curve calculated with help of (|33[) . Therefore, it seems that this issue is not 
significant for a one Jupiter mass planet. It is even less important for more 
massive planets whose evolution time scales show a stronger dependence on 
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Fig. 5 a) The evolution time scale t ev as a function of Pi in units of 1 day. The planet's 
mass M = lMj. Curves are numbered below in order of increasing magnitude at a fixed 
value of Pi for Pi > 2. The solid curve (3) corresponds to t ev obtained with help of H56H 
together with AE calculated for a non-rotating star. The short dashed curve (1) is for the 
retrograde case with fl(t) given by H49I I. The corresponding plot for the prograde case almost 
coincides with the solid curve, and so is not shown. The long dashed (4), dot dashed (5) and 
dotted curves (6) are for tides excited in the star (AE p [ = 0). They are for the retrograde, 
a non-rotating and prograde cases, respectively. The dot dot dashed curve (2) is for the 
retrograde case with all tides , but with the energy transfer AE reduced by a factor of 1.5 as 
compared to the value given by H33I I. b) As in a) but for M = 5Mj and with the difference 
that in this case the long dashed (5) , dotted (4) and dot dashed curves (6) determined from 
tides exerted in the star alone correspond to the retrograde, non-rotating and prograde 
cases respectively, c) The ratios of the evolution time scales corresponding to the prograde 
and retrograde cases, tf„ and t r ev , to the time scale corresponding to the non-rotating star, 
t® v , are plotted as functions of Pi. The solid and dashed curves respectively apply when 
tides exerted in the star and planet are both included. The dotted and dot dashed curves 
respectively apply when only tides in the star are included. 
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the relative direction of the orbital motion and of the stellar rotation, see e.g. 



Fig. 5(b) 



6.2.2 Conditions leading to stochastic instability or effective dissipation of the 
mode energy 

As discussed in Kochanek (1992), Kosovichev & Novikov (1992), Mardling 
(1995)a,b, Mardling & Aarseth 2001, Ivanov & Papaloizou (2004), (2007), 
when the orbital semi-major axis is larger than some critical value, a st , stellar 
and planetary perturbations excited at successive periastron passages have 
approximately uncorrelated phases. This happens when the phase change of 
perturbed quantities associated with a normal mode occurring during one 
orbital period, that arises on account of the change induced in the orbital 
period itself, is larger than some critical number of order of unity. This change 
in the orbital period results from the energy interchange between the orbit 
and the normal modes. This number was estimated numerically by Ivanov & 
Papaloizou (2004) in the simplest case of only one mode interacting with the 
orbital motion due to tides. 

In the situation where the phases are uncorrelated the mode energy, which 
is proportional to square of the mode amplitude grows on average in proportion 
to the number of periastron passages. This effect is referred to as the stochastic 
instability. 

From the results obtained in Ivanov & Papaloizou (2004) it follows that 
the stochastic instability sets in when 

a>a st = {ue DT y^ *'* (5*) an, (57) 

where Co = uj/f2*, R p i is the planet's radius and Rj and Mj are the radius 
and mass of Jupiter, respectively. We use equation (I57[) to estimate a st , for 
definiteness setting Co = 0.5 there, which is appropriate for g modes excited 
in the star and m p i = Mj, R p i — Rj. Since both stellar and planetary tides 



depend on time t we plot the dependence of a st on t in Fig. 6(a) for different 
values of Pi = 1, 2, 3, 4 and 5. 

One can see from this figure that the critical semi-major axis a st is typically 
of the order 0.1 — lOcra. Note, however, that the analysis of Ivanov & Papaloizou 
2004 can somewhat overestimate the value of a st since it was done under the 
assumption that only one normal mode plays a significant role. For the problem 
on hand, several modes are excited in the star and and in the planet which 
could give a comparable contribution to the energy exchanged with the orbit. 
This is likely to increase 'the degree of stochasticity' in the system leading to 
smaller values of a s t- However, in any case it seems highly improbable that 
dynamic tides can operate in the stochastic regime at scales much smaller 
than, say, O.lau,, for the whole of the interesting range of Pi. Thus, provided 
that dissipation of mode energy is not unexpectedly large, dynamical tides are 
not efficient at such scales so that quasi-static tides need be invoked to explain 
further tidal circularisation. 
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b) 



Fig. 6 a) The scale a s t as a function of time t. The rotational history of the star is given 
by II49II . Solid curves correspond to the retrograde orbits with values of Pi = 1,2,3,4,5 
with curves having larger values at given moment of time corresponding to larger values 
of Pi. The difference between the prograde and retrograde orbits is not well pronounced. 
As an example the dashed curve shows the dependence of a s t on t for the prograde orbit 
with Pi = 5. b) Same as the upper panel but a s t is shown as a function of Pi for a time 
or stellar age of t = 10 7 yr. The solid and dashed curves show the retrograde and prograde 
cases, respectively. 
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p 



b) 



Fig. 7 a) The characteristic time of non-linear dissipation, t n l as a function of the dimen- 
sionless orbital period after circularisation, Pi for a star with rotational period 1 day. The 
solid and dashed curves correspond to the prograde and retrograde cases, respectively, b) 
Ratio of the energy transfer AE to the threshold energy E t ^ r for the parametric instability 
to operate. The value of E t ^ r is taken to be 10 s5 ergs, see the text for discussion. 
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Fig. 8 The evolution time-scale t ev calculated under the assumption that only tides exerted 
in the star are effective and that the energy transfer is proportional to the attenuation factor 
given by H59II . Note that here we define t ev as a time-scale of evolution of the semi- major 
axis from its initial value equal taken to be lOau to a final value taken to be 0.5au. The 
solid, dashed and dotted curves correspond to retrograde, non-rotating and prograde cases, 
respectively. The planet's mass M = lMj. 

However, in addition to the stochastic instability there is another poten- 
tially important process, which may facilitate the orbital evolution at large 
values of a. This is through the so-called non-linear dissipation of the energy 
E nm stored in the perturbed g-modes of the star (Kumar & Goodman 1996, 
hereafter KG). This results from the parametric excitation of low frequency 
high order normal modes through their non-linear interaction with the modes 
excited by tides. Provided that this energy is larger than a certain threshold 
value, E t hr, the characteristic time scale of dissipation of the energy, t n i is 
given by the expression (KG) 

//10 38 eros\ , s 
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When this time is smaller than than the orbital period P or h the energy trans- 
ferred during a single periastron flyby to the normal model can be dissipated 
in one orbital period, and, accordingly, the expression for AE can be simply 
added to describe the effect of successive periastron passages. On the other 
hand, when / = P or b/tni *C 1, assuming there is no build up, the energy 
transferred from the orbit to the star that is actually dissipated should be 
attenuated by a factor of the order of /. A more accurate calculation, which 
takes into account a possibility of resonances between the orbit and a mode, 
and therefore a build up transferred energy in the star, gives what we describe 
as an attenuation factor, f a , in the form (see Appendix C) 

fa = (l-cxp(-2/))(l+exp(-2/)-2cxp(-/)cos(27r W „/l2 orb ))- 1 , (59) 

which is valid for < / < oo. Here tu n is a typical value of eigen frequency of 
a mode excited by tides. For our problem it is of the order of 0.5 J?*. Note that 
f a — > 1 in the limit / — > oo as expected. In order to allow for the weakening of 
the orbital energy transfer in the limit / — > we adopt f a AE for the energy 
exchanged per pericentre passage in equation (|54p instead of AE itself. 



We show the dependence of t n i on Pi in Fig. 7(a) for a star with rotational 



period 1 day and the evolution time scales calculated taking into account the 



attenuation factor (|59[1 are shown in Fig. [5] As seen from Fig. 7(a) there is a 
marked difference between the cases of prograde and retrograde motion with 
the retrograde case corresponding to much smaller values of t n i at given value 
of Pi . This can increase the difference between probabilities of tidal capture of 
planets moving in opposite directions with respect to the stellar rotation. It is 
also important to note that the threshold value of the energy , E t hr, for our case 
can be estimated to be Ethr ~ 10 35 ergs. for a star with solar values of mass 
and radius and rotational period of order of a few daytd At sufficiently large 
values of Pi (see below) the amount of energy transferred during one periastron 
flyby of a planet moving on a prograde orbit is significantly smaller than , Ethr, 
making the non-linear dissipation of the mode energy impossible. This can 
further amplify the difference between the capture probabilities. Interestingly, 
the values of t n i for retrograde orbits are of the order of the orbital periods of 
planets with semimajor axes of order 10au - the scale, where the processes of 
tidal capture are expected to operate. 

In order to calculate t eV: as shown in Fig. [5J for simplicity we set AE p i = 
in (j5"4")) . One can see from this figure that when the attenuation effect is taken 
into account, dynamic tides can explain the initial stages of tidal circularisation 
up to Pi < 3 in contrast to the simpler case considered in the previous section, 
where the beginning of effective circularisation is expected up to Pi < 4, see 



Fig. 5(a) We also checked that typical values of AE used to calculate t ev 
are approximately equal to E t hr when Pi ~ 3. Thus, at larger values of Pi, 
non-linear dissipation of mode energy is likely to be unimportant. 



8 In this case KG estimate the value of E^r to be approximately 10 48 (a;t t h) — ^l^ergs, 
where u> is a typical eigen frequency of a mode undergoing parametric instability and t t h 
is the stellar thermal time. Taking u) ps 0.5J7* ~ 3 ■ 10 — 4 s~ 1 and tth ~ 3 X 10 7 j/r we get 
Ethr ~ W 35 ergs. 
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The results discussed in this section confirm our previously made asser- 
tion that dynamic tides can lead to an efficient decrease of semimajor axis on 
scales corresponding to tidal capture both as a result of the stochastic insta- 
bility and the possibility of non-linear dissipation of energy stored in normal 
modes excited in the star. However, at smaller scales, the evolution of a due to 
the processes we have discussed clearly stalls, necessitating the consideration 
of different nonlinear effects, such as wave breaking (e.g. Barker & Ogilvie 
(2010), Barker (2011)), operating on normal modes, or other phenomena such 
as turbulent friction operating on quasi-static tides to be able to describe later 
stages of orbital evolution. 



7 Discussion 

In this paper we have developed an analytic theory for calculating the tidal 
response of a slowly rotating star resulting from a fly by of a planetary per- 
turber in a parabolic orbit, with orbital plane at an arbitrary inclination to 
the stellar equatorial plane. Expressions for the energy transferred to the star 
were presented in section [3] (see equations (|2T|) -([55 | '). In order to develop these 
expressions only the first order corrections to the normal mode eigenfrequen- 
cies appropriate to the non rotating star were incorporated. This implicitly 
assumes that these modes only change slightly and continue to dominate the 
tidal response in this form. In order to confirm these results numerical solutions 
of the linear tidal problem were considered in 2] In these calculations Coriolis 
forces, but not centrifugal forces, were taken into account. Thus the equilib- 
rium model was spherical. Beyond this no assumption of slow rotation was 
made. In practice the calculations were done for a poly trope of index n = 3. 
Detailed results for prograde encounters, retrograde encounters and encoun- 
ters with orbital plane perpendicular to the equatorial plane were presented 
in section [5] These were compared with results obtained from the analytic 
treatment in section [51 In general the agreement between the two approaches 
was good. Both agreed that the tidal interaction was significantly more effec- 
tive for retrograde encounters. The energy exchanged between the orbit and 
the star was about one order of magnitude larger for retrograde rotation with 
fl = 0.1 ]?» compared to prograde rotation with the same angular velocity 
when r) = 4v2 ; with the deviation becoming larger as 7] increased. Nonethe- 
less the numerically obtained results showed an increasing deviation from the 
analytically determined ones for increasing angular velocity of retrograde ro- 
tation. The maximum departure was about a factor of two for = 0.25 and 
r] = 4\/2 increasing to a factor of 2.5 when 77 was increased to 8. However, 
due to the very rapid decline in the strength of tidal effects with increasing 
pcriccntrc distance, we note that such deviations are scarcely significant for 
practical applications. 

The time scale for the initial stage of tidal circularisation inferred from our 
results for an orbit with initial semi-major axis of lOait was given in section RT21 
for planets of lMj and 5Mj. The rotation of the star, taken to be one solar 
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mass, was assumed to follow the Skumanich (1972) law. When tides from 
both the planet were included, initial circularisation down to circular orbits of 
periods < 5days could be readily accounted for both masses. When the effects 
of tides in the planet were neglected, our results indicated that tides acting in 
the star alone could be capable of producing a significant change in the planet's 
semi- major axis for final periods < Adays for both masses. Provided that 
some additional processes such as quasi-static tides are sufficiently efficient at 
smaller scales our time scales may give an estimate of a characteristic time of 
the whole stage of circularisation. 



7.1 Application to the circularisation of Hot Jupiters 

The discovery of Hot Jupiters in orbits highly inclined to the stellar equatorial 
plane with some orbits even being retrograde has given support to the idea 
that at least some of these objects were scattered into highly elongated orbits 
that were the circularised by tidal effects. The process of the excitation of os- 
cillations by tides acting on the planets and the consequent circularisation has 
been discussed by Ivanov & Papaloizou (2004, 2007, 2010) and Papaloizou& 
Ivanov (2010). Here this work has been extended to consider tides excited in 
the star. These may be significant because in contrast to tides acting on the 
planet, they are not equally effective for prograde and retrograde rotation. This 
feature should be borne in mind when considering the distribution of objects as 
a function of inclination. Hebrard et al. (2010) point out that of the 37 systems 
for which the inclination between the orbital and stellar equatorial planes has 
been measured at the time of writing, about one third are strongly misaligned 
with the distribution of inclinations being dependent on the planet mass such 
that retrograde systems only appear for M <^ 2.5Mj. Our results indicate 
that under the assumptions we made about the stellar angular velocity, circu- 
larisation times for prograde and retrograde orbits should differ by a factor of 
two at least initially, for a 'typical' period after circularisation of the order of 
four days, and that systems should be seen up to the same limiting period in 
both types of system. There would thus be some asymmetry produced by the 
effect of tides in this case though this claim should be confirmed by numerical 
N-body experiments incorporating planet-planet scattering together with tidal 
effects that take account of effects resulting from stellar rotation. The asym- 
metry may be more pronounced for masses M >^ 5M,/. However, at present 
there are relatively few objects with measured orbital inclinations with respect 
to the stellar rotation axis measured in this mass range. Future observations 
may clarify this issue. Similarly results may differ for different histories of the 
stellar rotation frequency and/or different stellar masses. Winn et all (2010) 
point out that there is an increase in the fraction of planets on inclined orbits 
for stars with larger masses. These may have been rotating more rapidly in 
the past, thus increasing the asymmetry between the prograde and retrograde 
tidal interactions. 
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It is clear that the result for the energy transfer obtained in this paper 
can also be directly applied to a general situation where two rotating stars of 
comparable mass with misaligned directions of rotational axes tidally interact 
with each other. Therefore, our results may be of interest for the problem of 
tidal capture of stars in young stellar clusters. 

Our results can also be generalised in several ways. Firstly, expressions for 
the transfer of angular momentum can be obtained by a calculation similar 
to what is made in this paper for the energy transfer. Note that all three 
components of the angular momentum ( with respect to the coordinate sys- 
tem defined with respect to the stellar rotation axis) are transferred in the 
case of a general inclined encounter. Secondly, a more realistic stellar model 
should be used to calculate the eigenfrequencies, overlap integrals and the fre- 
quency splitting coefficients (3 r . Thirdly, in order to reach better agreement 
between the theory and the numerical calculations, a better theory of the 
tidal excitation of high order (/-modes in a rotating star should, in principle, 
be developed. In the framework of such a theory one may try to relax the 
assumption of smallness of the stellar angular velocity in comparison to the 
mode eigenfrequencies. These issues will be addressed in future investigations. 
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Appendix A 

Explicit expressions for df£p n 

We list below the explicit expressions for the Wigner d— functions, which arc used in our 
calculations. We need 




(1 ± cos/3) 2 



4 




sin/3(l ± cosff) (2 ) 




sin 2 /3, (60) 



and the functions obtained from (160 H by the rule d\ 
Additionally, we use 



.(2) 

-n,±2 ~ 



for n = 0, 1,2. 




3 cos 2 p-1 
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= T 




— sin p cos p. 



(61) 
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Appendix B 

Transfer of angular momentum during the periastron fiyby 

The rate of angular momentum transfer, j, as well as its value of the periastron flyby, 
A J = J °° dtj are determined by the expression 

dV(rxVlT), (62) 

where r is the radius vector, and p is the density perturbation, which can be found from 
the continuity equation 

P = -V-(p£) = -V ^ (-^-(pr 2 &)-6pCs)Y 2 , n , (63) 
^— ' r r dr 

n 

where we use equations (15), (16) and the known properties of the spherical functions. Now 
we substitute equations (3) and 1 1631 ) in H62I I, integrate by part the term proportional to 
^r(P r2 £r) and use (20) to obtain 

j = Q^^nb fc Jfc,„, Jfc, n = / dnY 2ik (e r X rVY 2 , n ), (64) 



/i . k 



where integration is performed over the solid angle, and e r in the unit vector in the radial 
direction. The quantity 3 k n can be easily evaluated with help of the theory of the vector 
spherical functions (see e.g. Varshalovich et al 1989) with the result 



Jfc.n = i(-l) fc [n5 fcj _ n e z ( v /(2 - fc)(3 + k)6 k _ (n+1) e_i - ^(2 + fc)(3 - fc)5 M _„ei)], 

(65) 

where e±i = ^i(ej, ± ie,,), and e x , e y , e z are the usual unit vectors associated with the 
Cartesian coordinates (x,y,z). 

In order to calculate A3 we express the quantities b k and A„ in (1 6411 through their 
respective Fourier transforms and integrate the result over time using (7) to get 

AJ = Qj2 M n,kJk,n, M n<k = 2n(-l) n J dab k (a)(A- n (a))*. (66) 

n,k 

Now we substitute (21) in (1661) and calculate the integral over a in (1661) taking into account 
only the resonance contribution. The calculation is analogous to the one made for the energy 
transfer in Section 3 with the result 

M nM = i—Q[(-l) n A k (u+)A- n (u>+y - (-1)*!^+ )I_ fc (w+ )*]. (67) 

Now we substitute (1651) and (1676 in (1661) and get the angular momentum transfer in the form 
A3 = AJ z e z +AJ- 1 e- 1 +AJ 1 e 1 , (68) 

where 

= ^ZT" 12 fc (l^(^fc)| 2 ~ l^-fc(^J| 2 )' ( 69 ) 
k 

ZU_1 = 12 V( 2 - *)(3 + fc)(i fc K + )4 +1 K + )* + A_ (fc+1) (o;+ fe )I_ fc (a;+ fe )*, 

(70) 
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and AJ 1 = — (AJ" 1 )* . Note that to obtain the latter relation we make a change k — > —k 
in the summation series Jtjfc> P when considering the term proportional to ei . The Cartesian 
components of the vector of the angular momentum transfer, AJ X and AJ V , can be found 
from the relations AJ X = -^(AJ- 1 - AJ 1 ) and AJV = -^(AJ' 1 + AJ 1 ). We have 

AJX = ^IT~ 51 A/(2-fe)(3 + fc)(fl fc (w+)fl fc+ i(o;+) + J fc (w+)J Jk+1 (w+) 
k 

+ fl-fc(a>t fc )fl_( fc+ i)(w+ fe )+/- fc (a;+ jfc )/_ (fe+1) (a>+ fc )), (71) 
and ^ 

4 JH = 0[5L ^/(2 - fc)(3 + *)(J fc+ i(«+)flfc(w+) - 4(u>+)it fc+1 (u>+) 

+ /_( fc+1 )(^J J R- fc (^+ fc )-/_ fc (^+ fc ) J R_( fe+1) (^+ fc )), (72) 

where -Rfe(cr) and Ik(&) are the real and imaginary part of Afc(cr), respectively, and we 
remind that Afe(cr) is defined in equation (14). Equations 1691 1. 1171 I I and 1721 allows us to 
determine all three components of the angular momentum transfer after a tidal encounter. 



Appendix C 

Energy transfer from a highly eccentric orbit associated with a nor- 
mal mode with general damping rate 

Here we consider the energy exchange associated with an orbit of high eccentricity where 
the mode damping rate may be small. In this case the orbit is assumed to be quasi-periodic 
and so the system undergoes multiple close approaches in contrast to the single encounter 
considered in parabolic limit that applies when the mode decay time is shorter than the 
orbital period. The orbital period is assumed short enough so as to avoid the stochastic 
regime (see Ivanov & Papaloizou 2004). We obtain an estimate for the factor by which 
results for the parabolic case need to be multiplied in order to apply to this case. 

We begin by writing down equation 1181 governing the response of a single mode am- 
plitude allowing for mode damping (see Ivanov & Papaloizou 2004) in the form 

K +u>o6„ +2( 7n +iuj 1 )b n +2in-y n ^b n = S, with S = J d z xp{i* n -\/U), (73) 

where f n is the mode damping rate. Noting that the problem is very similar to that of a 
forced, damped harmonic oscillator, we identify the mode energy as 

£n = i(|M 2 +o; 2 |M 2 ). (74) 
Its rate of dissipation through mode damping follows from 1731 as 

— f = -2 7 „|fe n | 2 - in ln Q{b n (b n )* - (6„)*6n), (75) 
dt 

see e.g. Ivanov & Papaloizou 2004. 

We consider forcing in a periodic orbit with period P or t = 2ir/f} or i,. The potential U 
will then have a discrete Fourier series representation in the form 



oo 

U = y Uk exp (— ifci? or ;,t), 



(76) 
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where the Fourier amplitudes are given by 



l-n/0 orb 

U k = -^- I Uexp(ikn orb t)dt. (77) 

-7v/n orb 



2- 



Note that as the period tends to oo, or Q orb — > 0, with kQ orb remaining finite, 2-nlA k lfl orb 
approaches the continuous Fourier transform associated with a parabolic orbit. 

Because we consider an essentially linear response, equations II 731 -1 1 74 I I may be applied 
to each Fourier component separately and then the results combined. Denoting by S k , the 
quantity S evaluated for the corresponding term in the Fourier expansion for U, and the 
associated mode response amplitude by b nk , we have 

Kk = —, D = Jq - k 2 !7 2 rb + 2u^kn orb - 2i~{ n {kn orb - nO). (78) 

As we discuss below only contributions in summation series over k, where ujq ~ kft orb are 
important for our purposes. This means that the term proportional to 7 n J? in the expression 
for D in 1781 is much smaller than other terms as long as Q <C ujq. It is neglected below. 
In the same approximation the term proportional to fl in 1751 1 can be neglected as well (in 
fact these terms could be retained in the analysis but little is gained). 

Note that in order to get 17811 we assume that the mode damping rate is approximately 
constant. This is not formally the case when the dissipation is determined by non-linear 
processes considered in the text and -f n is a function of mode amplitude. However, from 
equation 17511 it follows that the characteristic time of mode amplitude decay is long com- 
pared with typical periods l/(kfi orb ) when -f n /(kfi orb ) <C 1 which corresponds to the case 
of interest. Then we can approximately consider 7„ to be constant and to correspond to 
the initial rate of of mode decay, 1 / t„i , when dealing with the Fourier harmonics with high 
value of k that have frequencies corresponding to those of the excited normal modes. 

We now evaluate the rate of energy dissipation by applying equation 1741 1 for each Fourier 
component and summing over k to give 

ds n -A k 2 n 2 ,\s k \ 2 

dt ln [[(fcfi ort - cl -Wo)^! - W0 -fefi orb ) + (a ; l)2]2 + 4 7 2fe2fi2 ]' V > 

k — — oo 

We may make an estimate for the above sum by noting that for a weakly damped oscillation 
mode, most of the contribution to the sum is expected to come from values of k corresponding 
to near resonance between a harmonic of the orbital period and a mode oscillation period. 
Neglecting the inconsequential quantity (o;^) 2 , which has been assumed to be negligible, 
resonance occurs when k is such that kQ orb = oj^ ± ujq, which in fact corresponds to two 
modes, one with positive frequency and the other with negative frequency. As the important 
values of k differ in sign, these can be dealt with separately and the results summed. Without 
loss of generality, we shall restrict attention to the positive frequency mode for which values 
of k near lo / Q orb S> 1 are expected to dominate the sum (we recall that u> = u)q + ) 
. Assuming other quantities vary smoothly with k, we replace k by fco, being the closest 
integer to ui/Q orb , everywhere except in the first bracket in the denominator of 1791 1 to 
obtain 

d£n lnU 2 \S ho \ 2 y> 1 

k — — oo 

The sum over k is readily performed by standard complex variable techniques to give 

d^ = _^ 2 \s k0 \ 2 ^ (gi) 

dt 2ujQ^2 orb 

where 

jr _ (exp(47T7 n /J7 or6 ) - 1) 

[exp(47T7 n /J7 ori) ) + 1 - 2cxp(27r7 n /J? orb )cos(27ra;/J? ori) )] 
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To interpret (1816 we remark first that when 27T7 n / J? or j, is large corresponding to a damping 
time shorter than the orbital period, T = 1. Then we can integrate H81II over an orbital period 
to obtain A(£ n ) oc \Sk \ 2 / ^ rb , which ultimately depends only on the continuous Fourier 
transform of the perturbing potential. This limit corresponds to a parabolic encounter. More 
generally T gives an estimate of the factor by which the energy exchange is modified when 
the damping time is longer than an orbital period and the motion can be regarded a quasi 
periodic. When 2ir r ) n / f! or ), is small, we have 

t ~ ^;l Q r "\ . ( 83 ) 

sin (7r£j/17 or (,J 

This factor could be large very close to resonances where u = kQ or i,, for some integer k. 
However, between resonances T has a minimum value •y„P or t>/2, while a reciprocal mean 
over mode frequencies would suggest a value r y n P rb- As orbital evolution very close to 
resonances is expected to be rapid, the latter values are expected to be more appropriate for 
estimating evolution times. Recall that for comparison with section [6.2.21 the mode decay 
rate "f n adopted here here should be identified with used there. 
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